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无 人 机 影像 三 维 重建 在 沙丘 形态 监测 中 的 应 用 
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摘 要 : 沙丘 形态 监测 对 开展 土地 沙化 防治 工作 其 有 重大 意义 。 应 用 无 人 机 ,分 别 在 2016 年 7 月 、 
2017 年 3 月 7 月 和 9 月 对 青海 湖 克 土 沙 区 4 个 沙丘 进行 影像 采集 ,借助 PhotoScan 软件 进行 三 维 
重建 生成 高 分 辩 率 的 沙丘 正 射 影像 和 DEM 数据 ,使 用 ArcGIS 软件 获取 沙丘 形态 参数 并 分 析 变 化 。 
实验 结果 表明 :基于 无 人 机 影像 监测 沙丘 形态 ,具有 可 行 性 和 准确 性 , 正 射影 像 精 度 较 高 ,DEM 数 
据 通过 校正 后 可 正确 反映 沙丘 形态 ,为 逐 天 高 精度 监测 变化 提供 了 可 能 。 
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沙丘 是 风沙 地 貌 学 的 主要 研究 对 象 ,关于 沙丘 
形态 的 形成 发 育 以 及 分 布 规律 是 其 重要 的 研究 内 容 
之 一 ，。 而 要 对 这 些 内 容 进行 研究 ,首先 需要 对 地 
貌 形态 进行 监测 ,获取 形态 数据 "“” 。 沙 丘 形 态 监 
测 研究 的 进展 ,大 致 可 分 为 4 个 阶段 "7, 监测 方 
法 ,也 由 最 开始 的 定性 描述 方法 ,发 展 为 利用 插 
标杆 法 等 的 定位 半 定 位 观测 方法 ” ,借助 卫星 和 雷 
达 影 像 的 遥感 监测 方法 “ ,以 及 近期 被 广泛 使 用 
的 利用 GPS .全 站 仪 和 三 维 激 光 扫 描 仪 等 的 三 维 测 
量 方法 ” 。 方 法 的 逐步 发 展 , 都 是 源 于 对 更 高 精 
度 沙 丘 形态 数据 的 追求 。2011 年 以 后 ,由 于 民用 及 
消费 级 无 人 机 的 普及 ,无 人 机 遥感 (unmanned aerial 
vehicle remote sensing，UAV-RS ) 技术 被 广泛 应 
FS! 。 相 较 于 传统 遥感 技术 ,其 具有 高 分 辩 率 、 
高 时 效 性 云层 下 成 像 和 移动 性 能 高 的 优点 。 近 年 
来 逐渐 开始 应 用 于 灾害 预报 监测 与 评估 , 正 射影 
BER ,及 三 维 重建 等 方面 "” 。 应 用 无 人 机 影 
像 进行 三 维 重 建 ,可 用 于 生成 数字 高 程 模型 DEM 和 
数字 地 表 模 型 DSMI] ,这 些 数据 可 用 于 文物 保 
护 四 1 ,灾害 评估 ,以 及 地 形 地 貌 等 方面 的 研 
究 -”*|。 应 用 无 人 机 影像 进行 三 维 重建 ,可 获得 高 
精度 的 DEM 数据 ,但 并 无 学 者 将 此 种 数据 用 于 
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沙丘 形态 监测 。 

本 文 应 用 大 性 公 司 的 精灵 4 无 人 机 在 青海 湖 克 
土 沙 区 对 典型 新 月 形 沙丘 进行 4 次 不 同时 期 野外 沙 
丘 影像 采集 ,并 在 室内 ,应 用 Agisoft 国人 公司 Photo- 
scan 软件 对 影像 进行 三 维 重 建 ” ,生成 正 射影 像 
和 数字 高 程 模型 DEM, 2454 Esi @ 公 司 的 数据 分 
析 软 件 ArcGIS™ 进行 沙丘 形态 参数 提取 与 比较 ,用 
于 监测 沙丘 形态 变化 。 尝 试用 这 种 方法 探究 应 用 无 
人 机 影像 对 沙丘 形态 进行 监测 的 可 行 性 和 准确 程 
度 。 从 而 实现 快速 准确 且 便 捷 地 获取 沙丘 形态 数据 
及 了 解 其 变化 情况 ,为 土地 沙化 防治 提供 支撑 。 


1 研究 区 概况 与 数据 获取 


1.1 研究 区 概况 

青海 湖 位 于 青藏 高 原 东北 部 , 北 依 大 通 山 , 南 傍 
青海 南山 , 东 靠 日 月 山 , 西 与 哈 拉 湖上 毗邻 ,构成 独立 
封闭 的 青海 湖 倪 地 地 貌 。 且 该 地 区 地 处 常年 西风 
华 、 东 部 季风 区 和 青藏 高 原 季 风 区 的 交汇 地 带 。 受 
地 形 海拔 和 大 气 环流 影响 ,形成 太阳 辐射 强烈 AT 
温 日 较 差 大 .区 域 风速 强劲 等 气候 特征 '” 。 年 平均 
气温 -0.2 ~0.6?C ,昼夜 温差 大 。 年 降水 量 380 ~ 
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420 mm ,年 蒸发 量 为 1 422.5 ~2 066.2 mm! , H 
特殊 的 地 理 位 置 ,造成 了 该 区 生态 环境 的 特殊 性 与 
脆弱 性 ,土地 沙漠 化 危害 更 加 严重 ,特别 是 东 岸 分 布 
有 大 片 的 沙丘 。 因 此 对 本 地 区 土地 沙化 问题 的 研究 
具有 维护 当地 生态 平衡 的 现实 意义 。 

用 于 沙丘 形态 监测 实验 的 4 个 沙丘 位 于 青海 湖 
东 岸 , 海 受 湾 东 南部 的 克 土 沙 区 ,地理 位 置 坐标 
为 36°%47'N,100°%47' 福 ,海拔 高 度 在 3 200 ~3 300 m, 
面积 约 0.62 km*( 图 1)。1 号 和 4 号 沙丘 为 典型 新 
月 形 沙丘 ,相对 高 度 分 别 为 36 和 30 m, WERE 
无 植被 覆盖 , 仅 在 丘 间 地 有 沙棘 和 沙 草 分布。2 号 
为 新 月 形 沙 礁 ,相对 高 度 45 m, 丘 间 地 有 沙棘 和 沙 
蒿 分 布 。3 号 为 大 型 沙 垄 ,相对 高 度 约 100 m, HR 
北 侧 存 在 小 幅度 波动 。 

1.2 技术 路 线 

本 文 应 用 无 人 机 对 沙丘 地 表 进 行 三 维 重 建 的 总 

体 技术 路 线 如 图 2 所 示 。 
1.3 无 人 机 数据 获取 方法 
通过 人 工 操控 大 疆 精 灵 4 无 人 机 的 方式 进行 影 
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图 1 研究 区 沙丘 位 置 无 人 机 影像 (2017 年 7 H) 
示意 图 (1 : 12 000) 

Fig.l Location of study dunes in unmanned 


aerial vehicle images , July ,2017 (1 : 12 000) 
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青海 湖 湖 东 沙 地 
4 个 样本 沙丘 


无 人 机 遥感 差分 GPS 沙丘 
ee 沙丘 数字 高 程 
像 ) 过 模型 DEM 


样本 个 体 长 : PATHE 
度 、 宽 度 高 人 线 位 置 及 长 
度 和 体积 A 度 
样本 个 体 的 形 
态 特征 及 变化 
情况 
图 2 技术 路 线 


Fig.2 Technical route 


像 采 集 , 相 机 为 DJI FC330 ,焦距 4 mm。 控 制 无 人 机 
飞行 高 度 为 120 m, 飞 行 速 度 在 3 mss 左右 ,影像 
拍摄 间隔 为 3 s ,拍摄 时 相机 镜头 垂直 向 下 。 尽 量 保 
UEFA TESA KO A [el Be BE A SF I] Be BE LE 70% 以 
上 ,以便 后 期 进行 地 表 三 维 重建 能 生成 更 精确 数据 ， 
理论 上 是 航向 重奏 60% , Hh] EA 40% ”” 。 本 
文中 ,在 2016 年 7 月 15 16 日 ,天 气 晴 转 多 云 ;2017 
年 3 月 24~26 日 ,有 风 , 风 速 2~6m.s ;2017 年 
7 月 23 ~25 日 , 睛 ;2017 年 9 月 13 ~15 日 ,上 晴 。 分 
别 对 4 个 沙丘 进行 实验 区 无 人 机 影像 采集 。 其 中 1 
号 2 号 和 4 号 沙丘 的 四 期 影像 全 部 成 功 获取 ,3 号 
沙丘 仅 成 功 获取 3 期 影像 ,2017 年 3 月 的 无 人 机 影 
像 因 采集 区 域 位 置 错误 而 无 法 使 用 。 
1.4 三 维 重建 

对 无 人 机 采集 影像 进行 三 维 重建 ,并 最 终生 成 
所 需 的 DEM 数据 和 正 射 影像 。 三 维 重建 的 基本 原 
SHE All FAW = fae ,已 知 参数 (图 3) 为 : 基 
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图 3 三 维 重建 的 基本 原理 


Fig.3 Foundation of 3D reconstruction 


线 (相机 两 次 拍摄 时 镜头 的 连 线 ) 长 ,相机 焦距 和 相 
机 CCD 大 小 ,目标 物体 C 在 两 个 相机 CCD 上 的 位 
置 。 利 用 相似 三 角形 可 以 建立 方程 , 求 出 目标 物体 
C 距 基 线 的 距离 以 及 其 距 相 机 两 次 拍摄 时 镜头 的 距 
离 。 实 际 中 要 解 算 出 目标 物体 C 在 空间 中 实际 位 
置 ,还 需要 相机 姿态 参数 及 相机 拍摄 时 的 GPS 位 置 
等 信息 。 

在 相 邻 两 次 拍摄 的 无 人 机 影像 中 ,会 有 许多 类 
似 于 目标 物体 C 一 样 的 点 。 借 助 Photoscan 软件 ,可 
以 解 算 出 所 有 这 些 点 的 位 置信 息 。 使 用 软件 时 , 首 
先 将 沙丘 影像 的 照片 导入 ,之 后 依次 选择 “对 齐 照 
片 ”“ 生 成 网 格 ” 和 “生成 纹理 ”对 影像 进行 处 理 。 
这 些 步 又 中 涉及 的 参数 均 按 默认 设置 。 完 成 后 , 即 
可 输出 正 射影 像 及 DEM 数据 ,设置 坐标 系 为 WGS 
1984 UTM zone 47N, 
15 沙丘 形态 数据 提取 

沙丘 形态 参数 包括 沙丘 兰 线 长 度 ,位置 ,沙丘 体 

让 ,沙丘 丘 顶 海拔 高 度 , 丘 顶 高 度 变 化 ,沙丘 蚀 积 变 

化 位 置 及 体积 。 沙 丘 刨 积 变 化 是 指 不 同时 期 ,沙丘 
不 同 区域 高 度 变化 的 差 值 ,进行 比较 的 区 域 是 固定 
的 。 结 果 包 括 蚀 积 变化 图 ,反应 蚀 积 变化 位 置 和 总 
体 蚀 积 变化 体积 。 

借助 ArcGIS 软件 ,对 正 射 影像 和 DEM 数据 进 
行 沙丘 形态 数据 提取 。 沙 丘 个 体 的 兰 线 长 度 ,位 置 
及 沙丘 边界 等 形态 参数 可 通过 对 正 射影 像 进行 矢量 
化 的 方式 获得 。 

沙丘 体积 , 则 应 用 3D Analyst Tools-Functional 
Surface-Surface Volume 工具 求 得 。 输 入 选择 沙丘 的 
DEM 数据 ,基准 面 高 度 选择 沙丘 DEM 数据 的 最 低 
值 ,或 丘 间 地 高 度 。 选 择 Above ,计算 基准 面 之 上 的 


体积 。 

之 后 ,应 用 栅 格 计算 器 Raster Calculator 对 两 次 
不 同时 期 同一 沙丘 的 DEM 数据 做 差 计算 ,从 而 得 知 
沙丘 形态 变化 情况 ,0 代表 不 发 生变 化 , 正 值 代 表 堆 
积 , 负 值 代表 侵蚀 。 由 于 沙丘 形态 变化 在 不 同 部 位 
堆积 或 侵蚀 不 一 样 。 借 助 Surface Volume 工具 , 设 
置 基 准 面 为 0, 即 不 发 生变 化 。 参 数 选 择 Above 即 
可 求 得 沙丘 变化 时 发 生 堆 积 部 分 的 体积 ,选择 below 
即 可 求 得 沙丘 变化 时 发 生 侵蚀 部 分 的 体积 。 


2 无 人 机 数据 精度 检验 


应 用 Photoscan 软件 ,对 无 人 机 影像 进行 三 维 重 
建 ,生成 的 正 射影 像 和 DEM 数据 结果 ( 表 1) ,其 中 
2016 年 7 月 数据 经 过 稀 蚊 点 云 加 密 处 理 ,精度 较 
高 ,后面 三 期 数据 由 于 数据 量 加 大 ,未 进行 稀 玻 点 云 
加 密 处 理 。 正 射影 像 和 DEM 数据 均 是 通过 三 维 重 
建生 成 的 模型 获得 ,模型 的 空间 位 置信 息 会 决定 正 
射影 像 和 DEM 的 空间 位 置信 息 。 对 于 沙丘 形态 监 
测 ,准确 的 空间 位 置信 息 十 分 重要 ,因此 在 下 面 对 生 
成 数据 进行 检验 。 
2.1 水 平 精度 

应 用 合 众 思 壮 (UnionStrong ® ) 公司 的 集 思 
G10A 星 站 差分 RTK 记录 沙丘 的 空间 位 置信 息 。 记 
录 方 式 选 择 水 平 误差 和 高 程 误差 为 厘米 级 的 固定 
解 . 采 点 间隔 约 3m。 路 线 选择 沙丘 的 边界 以 及 兰 


i 


表 1 无 人 机 影像 生成 数据 结果 
Tab.1 Reconstruction 's result from unmanned 
aerial vehicle images 
影像 采  ” 正 射 影像 DEM 


采集 时 间 DEO MEEO ARE DHR 
/ tk /m /m 
2016 £7 H 15 68 0.016 0.032 
25 52 0.027 0.053 
3 号 96 0. 035 0.069 
4 号 71 0. 008 0.016 
2017 年 3 月 15 177 0.023 0.504 
2 号 233 0.025 0.482 
3 号 一 一 一 
4 号 91 0.023 0.236 
2017 年 7 A 15 167 0.017 0.281 
2 号 158 0.011 0.141 
3 号 251 0.028 0.374 
4 号 56 0.016 0.258 
2017 年 9 月 15 190 0.020 0. 235 
25 418 0.020 0.270 
3 号 427 0. 038 0. 392 
4 号 105 0.024 0.269 
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线 ,同时 ,在 1 号 2 号 和 4 号 沙丘 还 选择 沿 迎风 坡 
和 背风 坡 采 集 ,采集 点 数量 见 表 3。 从 采集 的 GPS 
数据 中 ,选取 沙丘 峭 线 部 分 采集 点 。 利 用 ArcGIS 处 
理 无 人 机 DEM 数据 提取 沙丘 疹 线 。 对 比 点 线 距 离 ， 
得 到 无 人 机 影像 三 维 重 建 数 据 的 水 平 精度 检验 结果 
(图 4、 表 2)。 


表 2 沙丘 峭 线 位 置 水 平 精度 检验 表 
Tab.2 Accuracy test of dune’s ridge 


line’s horizontal coordinates 


沙丘 编号 2017 年 7 月 
1 号 GPS 采集 点 数 14 
点 线 距 离 平 均值 /m 0. 85 
2 号 GPS 采集 点 数 42 
点 线 距 离 平均 值 /mm 1.10 
3 号 GPS 采集 点 数 83 
点 线 距离 平均 值 /m 2.09 
4 号 GPS 采集 点 数 17 
点 线 距离 平均 值 /m 1.42 
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应 用 差分 GPS 数据 采集 的 2017 年 7 HWER 
线 位 置 数据 与 2017 年 7 月 基于 无 人 机 DEM 数据 提 
取 的 浓 线 较为 接近 (图 4)。 通 过 在 ArcGIS 中 计算 
得 到 的 点 线 距 离 的 平均 值 及 涉及 GPS 采集 点 数量 
如 表 2 所 示 , 点 线 距离 平均 值 均 小 于 2.1 m, 其 中 1、 
2 和 4 号 沙丘 疹 线 位 置 点 线 距 离 平 均值 小 于 1.5 m, 
而 3 号 沙丘 点 线 距 离 平 均值 为 2.09 m。 考 虑 到 在 
野外 利用 差分 GPS 采集 位 置 数据 时 的 人 为 误差 和 
差分 GPS 自身 误差 影响 ,可 以 初步 证 明 应 用 此 种 方 
法 生成 的 正 射影 像 和 DEM 数据 在 空间 水 平 坐 标 上 
精度 较 高 。 
2.2 BERE 

应 用 与 采集 无 人 机 影像 相同 时 间 采 集 的 差分 
GPS 数据 ,在 ArcGIS 软件 中 将 点 数据 转化 为 机 格 数 
据 , 栅 格 数据 格 网 大 小 与 无 人 机 影像 生成 DEM 数据 
一 致 。 上 一 小 节 已 论证 无 人 机 影像 三 维 重建 数据 在 
水 平 坐标 上 精度 较 高 ,因此 可 将 差分 GPS 经 处 理 得 
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图 4 水平 精度 检验 示意 图 (2017 年 7 H) 


Fig.4 Accuracy test of data’s horizontal coordinates , July ,2017 
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%3 DEM 高 程 数 据 与 差分 GPS 高 程 数据 差异 表 
Tab.3 Elevation data’s difference between 
DEM and differential GPS 


沙丘 2016 Æ 20174 20174 2017 年 
编号 7 月 3 月 7 月 9 月 
1 号 GPS 采集 点 数 116 731 798 625 
高 程 差 极 大 值 / m 116.7 101.7 107.2 69.6 
高 程 差 极 小 值 / m 90.9 97.9 101.0 42.4 
高 程 差 平 均值 / m 95.8 100.0 104.1 57 
高 程 差 方 差 12.7 0.6 1.9 9.0 
2 号 GPS 采集 点 数 44 423 579 502 
高 程 差 极 大 值 /m 92.2 109.5 141.9 67.4 
高 程 差 极 小 值 / m 86.3 94.0 100.8 59.7 
高 程 差 平 均值 / m 89 103.3 120.7 63.8 
高 程 差 方差 3.2 9.1 13.1 1.2 
3 号 GPS 采集 点 数 91 - 380 128 
高 程 差 极 大 值 / m 91.7 - 101.1 80.1 
高 程 差 极 小 值 / m 80.6 z 84.9 68.2 
高 程 差 平 均值 / m 87.4 - 95 74.9 
高 程 差 方差 9.8 - 14.6 6.9 
4 号 GPS 采集 点 数 65 279 380 313 
高 程 差 极 大 值 / m 98.6 110.1 106.8 76.9 
高 程 差 极 小 值 / m 87.1 105.8 101.2 70.1 
高 程 差 平均 值 / m 92.8 108.2 104.5 73.7 
高 程 差 方差 7.3 0.8 0.9 1.9 


到 的 栅 格 数据 与 无 人 机 影像 三 维 重 建 得 到 的 DEM 
数据 进行 高 程 上 的 对 比 。 对 比 二 者 在 高 程 上 的 差异 
( 表 3) ,进行 DEM 高 程 数据 精度 检验 。 

应 用 无 人 机 影像 进行 三 维 重 建生 成 的 DEM 数 
据 在 高 程 数据 上 与 差分 GPS 数据 存在 较 大 偏差 ( 表 
3)。DEM 数据 均 高 于 差分 CPS 数据 。 但 二 者 差 值 
均 集 中 于 某 一 区 间 内 , 且 高 程 差 方 差 较 小 ,数值 波动 
有 一 定 规律 性 。 因 此 可 应 用 ArcGIS 中 的 顶 格 计算 
器 工具 ,对 DEM 数据 进行 处 理 , 让 各 沙丘 DEM 数据 


分 别 减 去 其 相应 的 平均 差 值 。 消 除 高 程 差 对 精确 提 
取 沙 丘 形态 参数 的 影响 。 

无 人 机 影像 三 维 重 建生 成 的 DEM 数据 的 高 程 
信息 远 高 于 差分 GPS 数据 的 高 程 信息 主要 是 因为 
无 人 机 DEM 数据 高 程 信息 是 基于 拍摄 时 无 人 机 机 
身上 的 GPS 传感器 记录 的 坐标 进行 运算 生成 的 ,而 
所 使 用 的 大 疆 精 灵 4 无 人 机 ,其 GPS 传感器 在 高 程 
上 的 误差 较 大 , 且 并 未 购买 可 以 提高 无 人 机 GPS 精 
度 的 地 面 基 站 ,因此 高 程 上 出 现 误差 。 

综 上 ,基于 无 人 机 影像 生成 的 正 射影 像 和 DEM 
数据 ,在 水 平 坐标 上 ,误差 较 小 。 而 在 高 程 上 面 , 误 
差 较 大 ,应 借助 辅助 高 程 数据 进行 校正 。 


3 沙丘 形态 变化 分 析 


应 用 无 人 机 影像 进行 三 维 重建 生成 的 正 射影 像 
和 DEM 数据 ,经 过 校正 ,导入 ArcGIS 获取 形态 参数 
( 表 4, 图 5.6) ,进行 沙丘 形态 监测 。 以 1 号 沙丘 为 


例 。 
3.1 蚀 积 形态 变化 


应 用 无 人 机 影像 生成 的 高 分 辨 率 的 DEM 数据 ， 
可 以 获知 不 同时 期 沙丘 高 度 ,高 度 变 化 , 蚀 积 体积 变 
化 情况 ( 表 4, 图 5)。 图 5d 反映 的 是 2016 年 7 月 至 
2017 年 7 月 1 号 沙丘 的 形态 变化 情况 ,在 沙丘 的 西 
北 侧 , 背 风 坡 和 丘 间 地 交界 处 ,有 更 严重 的 侵蚀 ( 深 
蓝 色 ) 出 现 , 可 以 判断 是 产生 了 涡流 ,因此 造成 侵 
蚀 ,涡流 现象 也 出 现在 图 Sa 的 相同 区 域 。 这 也 证 明 
T DEM 数据 分 辨 率 较 高 , 且 较 为 可 人 和信。 而 在 图 Sb 
和 5d 中 ,位 于 沙丘 南 侧 ,有 一 块 较为 明显 的 侵蚀 区 ， 
这 一 区 域 由 于 地 处 1 号 沙丘 (图 5 中 所 示 沙 丘 ) 与 2 
号 沙丘 (1 号 沙丘 正 南 侧 ) 中 间 的 鞍 部 区 域 ,有 较 强 
的 峡谷 风 , 因 此 出 现 了 明显 的 侵蚀 区 。 


表 4 沙丘 形态 数据 
Tab.4 Attributive parameters of dune 

沙丘 编号 时 间 BARE ee 高 度 丘 顶 的 高 度 变 化 ”变化 体积 中 SN 

Am /m /m /m 堆积 体积 / m” 侵蚀 体积 /mm 
15 2016 年 7 H 43.64 402 111 3 269.01 - - = 

2017 年 3 H 51.50 325 143 3 263.45 -5.56 0.00 95 392.92 

2017 年 7 月 80. 66 215 519 3 264. 24 0.79 15 200. 30 5 088. 43 

2017 年 9 A 90. 59 263 320 3 264. 83 0.59 10 647.77 64 675.02 

2016 年 7 A 2017 年 7 月 一 一 一 一 4.77 0. 00 85 281. 06 

注 :沙丘 体积 蚀 积 变化 体积 指 前 后 两 次 采集 重 友 区 域 的 变化 情况 ,如 表格 中 2017 年 3 月 1 号 沙丘 堆积 体积 结果 ,为 2017 年 3 H DEM 数据 减 去 


2016 年 7 月 DEM 数据 后 经 处 理 后 得 到 
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图 5 沙丘 形态 变化 示意 图 (1 号 沙丘 , 底 图 为 2017 年 7 月 无 人 机 影像 ) 
Fig.5 Dynamic changing of dune morphology (example:dune NO. 1, base map :unmanned aerial vehicle images , July ,2017 ) 


由 于 区 土 沙 区 冬 春 季 主 要 为 西风 西南 风 和 西 
北 风 , 且 风 速 较 大 ,因此 在 图 5a 中 ,侵蚀 主要 集 
中 于 西 侧 ,而 在 东北 侧 侵蚀 量 较 少 。 图 Sb 所 反映 时 
期 以 西南 和 东南 风 为 主 ,但 由 于 风速 较 小 ,主要 以 堆 
积 为 主 , 集 中 于 沙丘 北 侧 。 在 夏季 ,以 西南 和 西北 风 
HE ,但 风速 较 小 ,因此 在 图 Se 中 ,侵蚀 主要 集中 于 
沙丘 西 侧 ,而 堆积 集中 于 沙丘 东北 侧 。 

总 体 而 言 (图 5d) ,1 号 沙丘 整体 发 生 侵蚀 , E 
要 受 冬 春季 较为 强劲 的 西风 和 西南 风 影响 。 侵 蚀 严 
重 区 域 集中 于 沙丘 西北 侧 (涡流 ) 和 南 侧 (峡谷 风 影 
响 ) ,总 体 侵蚀 量 见 表 4。 尽 管 在 2017 年 3 月 至 7 
月 有 少许 堆积 ,但 由 于 2016 年 7 月 至 2017 年 3 月 
这 段 时 期 侵蚀 量 较 大 。 因 此 总 体 上 发 生 侵蚀 , 且 在 
图 Sd 中 ,无 堆积 区 域 出 现 。 表 4 所 反映 的 1 号 沙丘 
1 年 蚀 积 变化 体积 ,为 2017 年 7 月 无 人 机 DEM 数 
据 与 2016 年 7 月 无 人 机 DEM 数据 对 比 得 到 结果 ， 
无 中 间 过 程 。 因 此 仅 有 侵蚀 体积 数据 。 


3.2 疹 线 变化 特征 

高 分 状 率 正 射 影像 ,可 用 于 闪 线 长 度 \ 位 置 及 其 
变化 等 沙丘 形态 参数 的 提取 ( 表 3, 图 6)。 总 体 来 
看 , 丘 顶 的 位 置 变 化 不 大 ,固定 在 几 条 痊 线 相交 的 地 
方 。 而 丘 顶 南北 两 侧 相 对 较为 平缓 的 疹 线 部 分 , 移 
动 频繁 , 北 侧 桨 线 在 2016 年 7 月 至 2017 年 7 月 向 
西 移动 ,之 后 向 东 移 动 到 9 H; ii pa MAE ZR HE 2016 
年 7 月 至 2017 年 3 月 向 东 移 动 ,之 后 改 为 向 西 ,但 
即便 到 了 9 月 ,其 位 置 也 并 未 退回 到 最 开始 2016 年 
7 月 的 位 置 。 从 养 线 来 看 沙丘 整体 由 东 向 东北 发 生 
偏转 ,上 且 向 西南 移动 。 

综 上 ,应 用 无 人 机 影像 对 1 号 沙丘 进行 形态 监 
测 生成 的 结果 ,经 过 校正 后 ,其 形态 变化 情况 反映 
出 :堆积 主要 集中 东北 至 东方 向 ,与 克 土 沙 区 输 沙 势 
和 输 沙 量 主要 集中 于 北 北 东 到 东方 向 的 情况 基本 一 
致 "9 。 这 说 明 利用 无 人 机 对 沙丘 地 表 进 行 三 维 重 
建 监 测 沙丘 形态 变化 具有 可 行 性 和 准确 性 。 
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图 6 沙丘 疹 线 变化 示意 图 (1 号 沙丘 ， 
底 图 为 2017 年 9 月 无 人 机 影像 ) 
Fig.6 Changing of dune crest line(example:dune No. 1 , base 


map:unmanned aerial vehicle images , September ,2017 ) 
4 讨论 


41 无 人 机 三 维 重建 结果 精度 
应 用 无 人 机 影像 进行 三 维 重建 ,对 沙丘 形态 进 
行 监测 。 数 据 的 精度 主要 取决 于 无 人 机 GPS 的 定 
位 精度 ,可 通过 设置 如 差分 GPS 原理 一 样 的 地 面 基 
站 的 方式 ,提高 其 定位 精度 。 当 然 , 也 可 以 通过 
直接 用 差分 GPS 在 实地 进行 测量 ,在 后 期 利用 Pho- 
toScan 软件 对 无 人 机 影像 处 理 的 第 一 步 “ 对 齐 照片 ” 
开始 前 ,加 载 这 些 测量 点 作为 地 面 控制 点 ,以便 
能 生成 更 高 精度 的 数据 。 本 文中 则 是 在 生成 数据 
后 ,利用 差分 GPS 数据 对 无 人 机 影像 三 维 重 建 结果 
进行 精度 检验 和 校正 。 借 助 差分 GPS 数据 对 无 人 
影像 三 维 重建 生成 的 DEM 数据 进行 校正 这 个 方 
法 是 可 行 的 ,有 助 于 提高 DEM 数据 精度 ,从 而 准确 
的 对 沙丘 形态 进行 观测 。 差 分 GPS 的 误差 可 通过 
背 助 水 准点 数据 进行 消除 。 
本 文中 所 设计 的 精度 检验 方式 ,也 可 以 在 今后 
的 实验 之 中 进行 完善 。 可 以 在 沙丘 上 随机 设置 9 个 
1 mx1 m 的 彩色 布 条 或 泡沫 板 , 之 后 再 用 无 人 机 进 
行 影像 采集 。 影 像 采集 完成 后 ,在 利用 差分 GPS 在 


彩色 布 条 或 泡沫 板 上 记录 位 置 数据 。 既 可 以 用 做 地 
面 控 制 点 ,辅助 三 维 重建 。 也 可 以 作为 精度 检验 的 
标准 ,检验 正 射影 像 的 精度 ,并 利用 软件 对 影像 进行 
校正 ,以 获取 更 高 精度 的 正 射 影像 和 DEM 数据 。 通 
过 对 比 利 用 差分 GPS 采集 的 准 线 高 度 位 置 数据 与 
无 人 机 影像 三 维 重 建 得 到 DEM 数据 提取 出 的 冰 线 
高 度 位置 数 据 对 比 , 检 验 正 射影 像 和 DEM 数据 的 精 
RE ,但 此 种 方法 需要 在 野外 利用 差分 CPS HEAR AFAR 
数据 时 ,有 和 较 高 的 准确 性 ,减少 人 工 采 集 误差 。 
4.2 沙丘 形态 变化 与 监测 

应 用 无 人 机 数据 来 监测 沙丘 形态 变化 , 相 较 于 
传统 遥感 手段 ,其 具有 很 多 优势 ,分辩 率 高 是 其 中 最 
主要 的 一 点 。 并 且 获 取 数 据 的 时 间 较 为 灵活 ,只 要 
野外 条 件 允 许 , 便 可 以 利用 无 人 机 采集 影像 。 不 同 
于 使 用 差分 GPS 测量 沙丘 形态 ,使 用 无 人 机 采集 数 
据 较 为 省 时 省 力 ,采集 0.02 km (2 号 沙丘 :200 m x 
100 m) 仅 需 半 小 时 , 且 通 过 运算 可 生成 分 辩 率 达 厘 
米 级 的 正 射 影像 和 小 于 1 m 的 DEM 数据 。 青 海 湖 
克 土 沙 区 ,在 夜间 风速 较 大 ,如 2017 年 3 月 25 日 至 
26 日 夜间 (18 :00 到 23 :00) ,小 时 平均 风速 在 6 ~ 13 
m* s ,最 大 时 达 13.2 m .s"!。 在 这 两 天 的 白天 ， 
我 们 均 对 1 号 沙丘 进行 无 人 机 影像 采集 ,分析 后 结 
果 表 明 , 在 迎风 破 ( 东北) ,最 大 堆积 高 度 达 1. 1 m, 
背风 坡 ( 西 ) 最 大 堆积 高 度 为 0.6 m。 这 为 对 同一 沙 
丘 逐 天 监测 提供 了 可 能 。 但 目前 的 是 结果 仅 能 达到 
厘米 级 ,对 较为 微观 的 沙丘 形态 变化 无 法 准确 监测 ， 
目前 只 能 应 用 差分 GPs” 。 

图 5a 明显 的 反映 了 涡流 现象 。 当 然 ,图 Sb 
和 5d 南 侧翼 角 受 峡 谷 风 影 响 发 生 较 大 侵蚀 还 需要 
结合 风速 风向 和 集 沙 量 数据 深入 分 析 。 但 实验 结 
已 经 能 正确 反映 出 沙丘 形态 变化 情况 ,说明 应 用 无 
人 机 影像 进行 三 维 重建 生成 的 正 射影 像 和 DEM 数 
据 , 经 过 校正 后 ,可 以 用 于 沙丘 形态 监测 , 且 较 为 准 
确 。 


5 结 


ole 
G> 


(1) 应 用 无 人 机 采集 青海 湖 东 岸 克 土 沙 区 影 
像 ,并 进行 三 维 重 建生 成 厘米 级 的 正 射影 像 和 分 辨 
率 小 于 1 m 的 DEM 数据 ,数据 精度 在 水 平 坐 标 上 与 
差分 GPS 记录 结果 接近 ,在 高 程 坐标 上 与 差分 GPS 
高 程 数 据 相差 较 大 ,但 可 以 借助 差分 GPS 高 程 数据 
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进行 校正 。 校 正 后 的 结果 能 正确 反映 沙丘 形态 。 

(2) 1 号 沙丘 在 2016 年 7 月 至 2017 年 9 月 这 
段 时 间 ,沙丘 整体 发 生 侵蚀 , 主要 集中 于 沙丘 西 侧 、 
南 侧 和 北 侧 与 丘 间 地 交界 处 。 丘 顶 的 位 置 变化 不 
大 ,位 于 4 条 不 同时 期 疹 线 相交 的 位 置 。 疹 线 由 东 
向 东北 转动 ,沙丘 整体 西南 方向 移动 。 

(3) 通过 应 用 对 无 人 机 影像 进行 三 维 重建 的 方 
法 ,可 生成 高 分 辨 率 正 射影 像 和 DEM 数据 , 且 该 方 
法 在 沙丘 形态 监测 的 应 用 方面 具有 一 定 的 可 行 性 和 
准确 性 。 为 逐 天 准确 监测 沙丘 形态 并 分 析 其 变化 提 
供 了 可 能 ,但 目前 无 法 应 用 于 较为 微观 的 沙丘 形态 
变化 监测 中 。 
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Application of UAV image based 3D reconstruction in morphological 


monitoring of dunes 


ZHANG Ming-yuan', ZHANG Deng-shan'”, WU Wang-yang', TIAN Li-hui>, ZHOU Xin? 
(1 Key Laboratory of Environmental Changes and Natural Disaster , Ministry of Education , School of Geography , 
Beijing 100875 ,China; 2 Qinghai Academy of Agricultural Forestry Sciences ,Xining 810016 , Qinghai , China; 
3 State Key Laboratory of Plateau Ecology and Agriculture , Qinghai University , Xining 810016 , Qinghai , China) 


Abstract: Morphological monitoring of dunes has an important value to control land desertification. Compared with 
using differential GPS to monitor dune’s morphology ,the method using small civilian drone was easier and had high- 
er reliability. This method could not only save time and labor but also produce high resolution results. And the ima- 
ges could turn out high resolution dune’s morphological data which was smaller than 1 m. However ,few studies mo- 
nitored dunes morphology by using unmanned aerial vehicle to collect images. In this paper ,unmanned aerial vehi- 
cle (UAV) images of four dunes in Ketu, eastern of Qinghai Lake, Qinghai Province, China were taken by DJI 
Phantom 4,in July 2016, March, July and September 2017, respectively. These images were 3D reconstructed by 
PhotoScan to produce high resolution digital elevation model (DEM) and orthophoto images. And also dunes’ mor- 
phological data was collected in the same time by using UniStrong G10A real time kinematic platform ( differential 
GPS) to test orthophoto images’ accuracy and corrected these images if there were some geometric error. For the 
DEM, differential GPS’s data was selected to test its accuracy. The results showed that these DEM’s elevation data 
were average 90 to 110 m higher than differential GPS’s data. But it could be fixed by using raster calculator tool in 
ArcGIS. Then through analyzing DEM and orthophoto images which were rectified through above steps by using Arc- 
GIS ,the dynamic change of dune’s elevation, length and location of crest line ,and dune’s volume could be under- 
stood intuitively. The results showed as follows: firstly ,the orthophoto image and DEM’s resolution were smaller than 
1 m. And the orthophoto images’ horizontal position was close to differential GPS’s data. But the DEM’s elevation 
was much higher than differential GPS’s altitude. This error could be revised and after that the DEM data could ac- 
curately reflect the true morphology of dunes. Secondly ,the dynamic changes of sample dune’s morphology charac- 
ters from July 2016 to September 2017 were consistent with previous studies. The erosion areas were mainly concen- 
trated in the west and south of the sample dune. It also concentrated in the border area of leeward slope and the area 
between sand dunes where it was in the north of the sample dune. Meanwhile ,the sample dune moved to the down- 
wind direction of main wind (southwest) and its crest line rotated from east to northeast. The position of the sample 
dune’s crest changed little. It was just at the junction of the sample dune’s four crest lines which were extracted 
from orthophoto images in different time , respectively. Thirdly , applying high resolution DEM and orthophoto images 
in morphological monitoring of dunes was feasible and accurate. It made daily monitoring possible. 


Key words: dune; morphological monitoring; UAV remote sensing; 3D reconstruction. 


